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Abstract 

We study the energy-momentum pseudo-tensor of gravitational wave, and examine the one 
introduced by Landau-Lifshitz for a general gravitational field and the effective one recently 



used in literature. In short wavelength limit after Brill-Hartle average, both lead to the same 

0" 1 



gauge invariant stress tensor of gravitational wave. For relic gravitational waves in the ex- 
panding universe, we examine two forms of pressure, p gw and V gw , and trace the origin of 

£T) ' their difference to a coupling between gravitational waves and the background matter. The 

> '■ 

. difference is shown to be negligibly small for most of cosmic expansion stages starting from 

00 ' 

inflation. We demonstrate that the wave equation is equivalent to the energy conservation 
equation using the pressure V„ w that includes the mentioned coupling. 

o 

(N ; PACS: 04.00.00, 04.30.-w, 04.30.Db, 04.30.Tv 

1. Introduction 

X 

H 

Relic gravitational waves (RGW) is believed to be generated during the inflationary stage of the 
expanding universe, and distributed as a stochastic background [TJ [21 EJ @]. Its amplitude, spectral 
index, running index, etc, are mainly determined by the initial condition, carrying a wealth of physical 
information of very early universe pp. As a possible channel for its detection, RGW can induce the curl 
type of polarization of cosmic microwave background radiation (CMB) |U [7J [U [9] . Since RGW is important 
for the early universe, one natural question is what are the energy density and pressure of RGW. Different 
definitions of the energy-momentum tensor of RGW lead to different equation of state for long wavelength 
modes, and thereby yield different predictions on its cosmological influences [EMU]. The field equation 
of RGW has been commonly employed as a solid ground to study the evolution of RGW in the expanding 
universe, and we shall examine which definition corresponds to the wave equation of RGW, providing 



* yzh@ustc.edu.cn 



1 



another perspective to the issue of energy and pressure of RGW. In most cases for physical systems, 
field equation and energy conservation equation are two different, yet equivalent, ways to describe the 
systems. As is known, the field equation can be derived from the equation of energy conservation, and 
vice versa, for such systems as oscillator, electromagnetic wave in vacuum, scalar field (both in flat and 
curved background spacetimes). The energy and pressure of these systems are well defined, and one can 
work with either equation for practical conveniences. For gravitational field, the problem of the energy- 
momentum is a thorny issue. As it stands, the energy-momentum of gravitational field can only be defined 
as a pseudotensor t a/3 , in stead of a tensor. This is due to its non-locality due to Einstein's equivalence 
principle [12J. The disadvantage of pseudotensor t a @ is that it is not unique. Various forms have been 
proposed for t a @ , each being based on certain specific considerations, such as the conservation of total 
energy momentum, symmetry, additivity, and so on (121 [151 [T4l [l5l [TBI [171 [151 [131 120] . 

In this paper, we shall study two definitions of the energy momentum pseudo-tensor (EMPT) of 
gravitational waves, emphasizing their relation to the wave equation. The paper is organized as follow. In 
section 2, we first examine Landau-Lifshitz' definition t a/3 of EMPT of a general gravitational field, and, 
for the simple case of for the flat spacetime, we demonstrate that the energy conservation of this t a @ leads 
to the field equation of gravitational waves (GW), and vice versa. In section 3, we investigate the short 
wavelength limit of Landau-Lifshitz' t a/3 , and apply Brill-Hartle average to it. In section 4, we examine 
the effective EMPT t^L of GW that appears in literature (20j EU [22l E]. In the short wavelength limit 
after a spatial average, t a f^ leads to the same expression as Landau-Lifshitz' does. In section 5, we study 
two definitions of the pressure, p gw and V gw , of RGW in the expanding universe, trace the origin of their 
difference and reveal its physical meaning and cosmological implications. Finally, in section 6, for RGW 
in the expanding Robertson-Walker (RW) spacetime, we examine the relation between the wave equation 
and the equation of energy conservation for the effective EMPT. 

2. Landau-Lifshitz' t a @ and Wave Equation 

For Landau-Lifshitz' definition of EMPT, we will explicitly show that the energy conservation is equiv- 
alent the field equation of GW in flat spacetime. For a general curved spacetime with the metric g a p, 
Landau-lifshitz' definition of the energy-momentum pseudo-tensor t a @ of gravitational field is given by [13J 

-G? Q Vg^ ,pg PM ,A + / Vg a " ,pg^ ,a) + 9 up 9x^ aX V ,P 
+I( 2 0«Yf - g^g^)(2g up g ar - g pa g VT )g VT , A g' CT ,„ (1) 

where g Q/3 = {—g) l ^ 2 g a ^, g = det \\ g a p ||= det \\ g a/3 ||, and the sub comma "," stands for the ordinary 
derivative. Einstein equation takes the form 

H a ^ v tltv = 16vr(- 5 )(r Q ' 3 + t aP ), (2) 
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where H a ^ u = g a ^g^ — g at, g^, antisymmetric in /3 and u, and T a/3 is the energy-momentum tensor 
of the matter. From Eq.([2]) it follows that 

[(-flXr^ + t°% = o (3) 

is identically satisfied. It implies that the four-momentum of the system, i.e, the gravitational field plus 
the matter, is conserved. This equation implicitly allows for interactions between the two components. 

For a general spacetime, the expression t a/3 in Eq.(TTJ) is rather complicated. We consider the flat 
spacetime background with the metric r]^ u = diag{— 1, 1, 1, 1) perturbed by GW, defined by h al3 = 
_g"/3 -j-^"/ 3 . Then the pseudo-tensor reduces to 

+i(2rf V M - v a ^V X n^Vu P Var ~ V P .Vur)h UT , x h pa ,„ (4) 

o 

As a physical quantity, GW has two degrees of freedom (polarizations), while the perturbation h a/3 has 
ten components, containing eight degrees of non-physical, gauge freedom. To remove these from h a/3 , we 
choose a synchronous coordinate system, in which 

h ? = 0, (5) 

and, furthermore, impose upon h 1 ^ the traceless-transverse (TT) condition as the following: 

h l t = 0, ft j j = 0. (6) 
Eq.© together with Eq.© are eight constraints. Then Eq.flU reduces to 

16vrG t af} = ^ V a/3 h ij >k h k itj + h a i>u h^' v - h a ijtfi'P - hP id h ij ' a 

+\h lha h 3 - j V a ^ x hi jiX . (7) 
In absence of matter T a " = 0, Eq.([3]) becomes (keeping up to the second order of h a ^) 

t aP ,p = 0, (8) 

i.e, the four-momentum of GW is conserved. We will check that this equation is consistent with the wave 
equation of GW. We take the derivative of t a(3 in Eq.©: 

-\h ijM hi^ - \h ij > X h^x- (9) 



For the "a = 0" component, it reads as 

= (io) 

From this equality it follows that the energy conservation 

t"% = (11) 



implies the wave equation of GW 



K ,„ = 0, (12) 



and vice versa. 

3. Short Wavelength Limit of Landau and Lifshitz' t a/3 

We are interested in GW in curved spacetimes, such as in an expanding universe, and the metric of 
spacetime is divided into two parts: 

9a/3 = 9$ +h Q p, (13) 

where g£l is the background spacetime metric, h a p is GW as small perturbations, 

I h a p |« 1. (14) 

As it stands, the general definition of t a ^ in Eq.([T]) is rather complex, and, in practice, it can be simplified 
for GW with wavelength being smaller than the scale of curvature of background spacetime. In short 
wavelength limit one can apply Brill-Hartle average [23l [12] 20]. 

There are three scales involved: the wavelength A of perturbations, the scale R of the curvature of 
background spacetime, and a scale L on which one takes average. The condition for the short wavelength 
limit is A <C L <C R. The background metric g^l is assumed to be a slowly varying function of spacetime, 
whereas the perturbation h a g changes significantly over small scales of space and time. Let A be the 
characteristic amplitude of h a p, A <C 1. In vacuum, the energy density of perturbations is ~ fj(4) 2, t ' ie 
curvature of background spacetime is ~ R~ 2 . If the curvature of background spacetime is induced by the 
energy density of perturbations, then, according to Einstein equation, one has R~ 2 ~ (§)(fy)(4) 2, ' e ' 
A ~ ^. By estimation of magnitude, one has 

o(°)~l ruL q {0) — 
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hfjiv ~ A, h^ Uta ~ A/X 



1 

R' 

By Eq. ([1]), t a ^ is quadratic in the first derivatives of the metric. The second order terms of perturbation 
in t a/3 , denoted as t^ a ^ , should contain such terms of the following form 

hh, hdh, dhdh, 

among them, the terms like dhdh are dominant, by estimation of the order of magnitude. So we will only 
keep the terms like dhdh in calculation. Then Landau-Lifshitz' t a " in Eq.([T]) is reduced to 

16ttG*°0 = A a P x A x \-A aX x A^ p + ^g Xp A x \A^ u 

-{g aX g^ v p A™ x + / P APil a) + 9 vp 9x»A aX V A^ p 

+\(2g aX g^ - g a Pg x »)(2g up g ar - g pa g VT )A VT \A^ p (15) 
o 

where 

a = {-g)- l/2 ^ , a = <T , a + ^ V CT <w (16) 
contains the derivatives. The terms of t a/3 in Eq. (jl5p are of the following form 

AA, ggAA, ggggAA. 

We expand A^ u a to second order of perturbation: 

J^iv ^ _ j^(0)fj,u ^ _|_ j^l)(iv ^ _|_ j^(2)^u ^ ^7) 

where 



A (0),v a = g Q»r a + 1^(0)^(0)^(0) ^ (18) 

A (l)^ a = _ h „u a + \g^gWp° hpa a - \gWr°gf) ja h^ - \g^ gf) t<x h»° , (19) 



A^ v a = \[g {0)p(T gtl a h^ + g^gfl^W" + g p %h^h^ 

-g^h^hp^ - g^ u h^h pc7 ^ + 2h^\ . (20) 

In short wavelength limit, we keep only the second order of perturbation, dhdh, in Eq . f|15|) . Only the first 
two terms of A^ 1 ' contain dh, which are denoted by 

^ a = _ h H» a + ±g(0)^g(0)fxr hpaa . (21) 

More specifically, dhdh actually come from a^a^. So we have 

WirGt^ = xa {1)Xp „ - a^ aX A a«^ „ + \g {0) ^ g^a^ P a^ „ 

-b ( ° )aA ^a (1) ^ ^ A + 9^9^- p a««« \\ 
+g(^g x y^ X u a^ p 

H2g^g^-g^g^)a^ x a^%. (22) 



Since the difference between the ordinary derivatives h^ u a and covariant derivatives with respect to 
background metric h pv \ a involves no terms like dh and will not contribute to the expected terms like 
dhdh, we replace the ordinary derivatives in Eq . f!22|) by covariant derivatives with respect to background 



(o) 

metric 5^, 



a^ w a = -h^ ]a + -g {0)fiv g {0)pcr h 



pa\a 



= -h^\ a + \g^h\ ai (23) 

where h = g^'^h^, and the subscript "1" denotes the covariant derivative using the background metric 
g£p, and the indices are raised and lowered with the background metric. In the gauge: h a ^ ^ = 0, 
h = h a a = 0, Eq^j reduces to 

a^ a = -h^ la , (24) 

and Eq. ([22~|) becomes 

+h Mu h« w + \h?°\ a h pa & - \g {0)al3 h^H pcTlfl . (25) 

This expression is still clumsy, and can be simplified in many applications using an averaging process 
adopted in Refs. [2UP23JP5Q]. Whenever regions of scale L are large enough to contain many wavelengths, 
it is natural and necessary to employ the average. The idea is similar to finding electric fields in a 
macroscopic dielectric. One is not interested in fine details of stochastic electric fields fluctuating over 
molecular scales, but only need to know electric fields averaged over a region large enough to contain many 
molecules and yet small compared with the size of the dielectric. In general, GW is similar to this. h pv 
consists of all sorts of modes of various wavelengths A, and those with A <C L are intractable, practically 
regarded as being stochastic. We will take average of Eq. (l25j) over a spacetime region of scale L. 

We only sketch the main idea of this kind of averaging procedure, omitting the technical details 
concerning the tensorial nature, which can be found in Ref.[20]. In right hand side of Eq. (i25|) there are 
six terms of the form dhdh. Consider a gradient H\ a composed of quadratic of dhdh. The average of 
H \ a over the region can be schematically defined as 



<H\ a {x) >= J H w (x')f(x,x')d 



J X , 



where f(x,x') is a weighting (window) function which falls smoothly to zero when \x — x'\ > L, and 
satisfies a normalization J f(x,x')d A x = 1. Then, 

< H\Jx) > 



]a (x) >= J[(Hf) w -Hf w )d A x'. 



The first term is actually a surface integral f s Hfd 3 x' which can be dropped as f\s = on the boundary 
surface, and the second term contains f\ a ~ f/L. Note that for GW, a derivative \ c ~ W" j\, so 



H\ c ~ if/A. Thus, the overall effect of averaging on the gradient Hi a is such that H/X — > H/L. Since 
we assume X/L -C 1, the divergence in average is vanishing < H\ a {x) >~ 0. From this result follow the 
rules of Brill-Hartle average (see §35.15 in Ref.[12j): 

• Gradients average out to zero; e.g., < (/i 7<5 \ a hiJ,v)\p >= 0; 

• One can freely integrate by parts, flipping derivatives from one h to the other; e.g., < h ryS h IMU \ aj3 >=< 
-/i 7<5 >; 

• Covariant derivatives commute; e.g., < h' yS h^ 1/ \ a p >=< K 1& h fa ,\p a >. 
The last rule is derived with the help of the equation of motion for h pu [20J 

h Hp e + 2R% p h<>° + R^K + RW h ° = o. (26) 

Applying these rule to Eq. (l25l) . only one term (l/2)h pa ^ a h pcr ^ remains, and we arrive at the energy 
momentum tensor for gravitational wave 

< >= t BH = J— < h^ a h pa V >, (27) 

327rG 

This expression has a desired property that it is a gauge invariant [5DJ [T2] and a physical observable. It 
has been commonly used in studies of gravitational radiation. For RGW in an expanding universe, Eq. (l27l) 
can be used, as long as the wavelengths under consideration are much shorter than the Hubble radius 
I/Hq ~ 3000Mpc. Otherwise, one need use Eq.([25|). 



4. The Effective t a / ff 

Another definition of the energy-momentum pseudo-tensor of GW has been used in literature, called 
the effective EMPT [2Ql Q21 EU E2l HTJ- We now examine its relevant properties. Consider a curved 
spacetime in Eq. ([T3l filled with matter. Einstein equation 

G a(3 = 8ttGT o13 (28) 

is expanded to the second order of perturbation, 

G (0)aP + Q(l)aP + Q(2)a/3 _ ^Q]^(0)aP + j.(l)a/3 + T (2)a/3] ^9) 

After rearranging some terms in the above equation we have 

= SirGlT^ + T& a P - -^-G^ af3 ] + [G^ - SttGT^}. (30) 

SttG 

The first order part of perturbation satisfies their field equation 

G (i)o0 = 8ttGT^, (31) 



which will be the main subject of cosmological perturbations. Then Eq. pop reduces to 

= 8ttG[T^ + T® a P + tf fJ \ , (32) 

where 

= -sb G<2> " s (33) 

is the effective EMPT [20J EE 122] ■ As the second order part, G (2)a/3 consists of quadratic terms of metric 
perturbations and will affect background spacetime, known as the "backreaction" effect [121 I2T1 1221 120j . 
Note that the second order matter term T^ a ^ also appears and affects the background spacetime as well. 
In practice, it is small and can be neglected. Explicitly, 

-8irGt a J f = -\g^g$R^ + \[g$h^ - g^h» v ]R^ 

+ h R (o)^ h ^ h ^ _ R$)h a Ph% (34) 

where the perturbed Ricci tensors are 

R (i)r> = g (o) mg (o)^ R {i) _ [g (o) mh »p + g (o)up hfia]R (o) (35) 

+ { R (0)na h vP + R (P)v* h tf] haj3 + R (°) h »c* h vP ( 36 ) 

with (see Eq. (35.58) in Ref.[T2]) 

Rfiv 2~( fy/if ^puv\ci ha(i\v "t" hav\fi )) (37) 

+h I/ a ^(h fia \i 3 — h^\ a ) 

-(h afi {p - hi\ a )(h Mv + h w]lx - Via)]- (38) 

Note that t"?f is not covariantly conserved, |/3 7^ 0, using the background metric. As it stands in 
Eq. (i34"l) . tgjj is different from the expression in Eq. (j25p by keeping only quadratic dhdh part of Landau- 
Lifshitz' t'fi, and is also clumsy for a general curved background spacetime. Similar to what we have 
done for Landau-Lifshitz' t a @ in the last section, one takes the average by applying the Brill-Hartle rule to 
Eq.flMD, yielding [20] 

<t ?/ >= 32b </lHa ^' /3> ' (39) 
where the gauge condition h a P |0 = 0, h a a = 0, are used. Thus, in the short wavelength limit after the 
Brill-Hartle average, the effective t°f^ defined in Eq. (f34]l is equivalent to Landau and Lifshitz' t a/3 , both 
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leading to the same expression. One sometimes extrapolates to use Eq. fl34p in the model of bouncing 
universe to investigate the backreaction of gravitons upon the background spacetime [11J. 

In the long wavelength limit [2T] 1221 E], the wavelength of perturbations is comparable to, or even 
larger than the scale of curvature of background spacetime. This can be relevant for study of the very 
early universe when the wavelengths of interest are even bigger than the horizon. Then the Brill-Hartle 
average procedure is not available. One may try to do spatial average of the EMPT in a fixed time slice of 
spacetime. However, there is a problem of gauge invariance with this method, and the issue is not settled 
down [SUES]. 

5. Pressure of RGW 



We consider a spatial flat homogeneous and isotropic RW spacetime filled with only RGW. The metric 

,(°) 

aft 



is given by Eq. (fT3l) with the background g^l = a 2 (r)diag(l,—l,—l,—l), and the perturbations in the 



synchronous coordinate are given by 



,0 o 

h a p= | | , i,j = 1,2,3. (40) 

—a 2 (T)hij 



In the following as RGW hij stands for what appear in Eq. (j40j) and satisfies the same TT condition 
h\ = dih l j = as in Eq.([6]). For this spacetime, a direct calculation of t°f^ given by Eq. (l34l) yields the 
energy density and pressure for RGW 

P 9W =4 = ^2lKh' kl h kl + l -{d m h kl d m h kl + h' kl h' k % (41) 

p gw = ~\i = ^G^il d mhkid m h kl - 5 -h' kl h' kl ], (42) 

where % = a'{T)/a{T) and the prime denotes the time derivative with respect to the conformal time. We 
drop the subscript " eff " for simplicity. In deriving the expressions of p gw and p gw , a total derivative has 
been dropped, respectively. No average is taken on the expressions (141!) and (l42j) , which is different from 
those given by the averaged expression in Eq. fl39p . As mentioned in last section, the effective t a ^ under 
discussion is not covariantly conserved. To get a conservation equation of RGW, Ref. [11] expands the 
Bianchi identity = to the the second order of perturbation: 

= = (G%.^ + (G"V ;M )« + (G%,^ 2 \ (43) 

requires each order should vanish respectively. Then the second order part 

(G m „;m) (2) =0, (44) 

which reads as 



dp g w 

dr 



+ m{ Pgw + V gw ) = 0, (45) 



where the pressure 

V ^^^m^^' (46) 

We will take a different approach to the conservation equation of RGW, and trace the origin of V gw . 
Taking covariant derivative on both sides of Einstein equation (f32|) using the background metric leads 
to 

G^ v ^ = ^G[T^ v +t^ w (47) 
By the Bianchi identity G^ ^ „l = 0, we have 

[T (f))li u + €]\ ft = 0, (48) 

which tells that the total energy and momentum tensor, i.e, the sum of matter and GW, is covariantly 
conserved, not tv itself. Only in absence of matter with T^ '^ v = 0, the EMPT of RGW is covariantly 
conserved with 

^ = °< ( 49 ) 
which, for the RW spacetime, explicitly reads as 



9pgw 



+ 3H{ Pgw +p gw ) = 0, (50) 



dr 

which is commonly used in literature on RGW. It is straightforward to show that the covariant conservation 
of the Brill-Hartle averaged tensor, 

t BH t\» = (51) 

also leads to Eq. (|50l) . We would refer to Eq. (l50|) as the "canonical" equation of energy conservation of 
RGW. It is noticed that Eq. ([45]) differs from Eq. (j50l) with V gw replacing p gw . 
In general, when the matter is present, T^°>^ v ^ 0, Eq. ([48j) gives 

^^ = - r( ° )M H^ (52) 

where the matter term — T^^u on the right hand side as an extra term can be interpreted as the 
interaction between GW and matter. We give a detailed calculation of this term. Starting with the 
covariant conservation equation of energy momentum tensor for matter, 

r<°>" „.„ = <), (53) 

expanding it to second order of perturbations, and taking spatial average, one has 

< „.„ >= T(°> v{fl + < r( 2 )^ ail > T^ a „- < T^ a Vii > rW" a = 0, (54) 

so that 

„!„ = -< ^ QM > T^ a v + < T^ a > T^ a . (55) 
10 



This result tells that T( q '''„| ]U represents a coupling between the metric perturbation and the background 
matter as given by the expression on the right hand side. Since the difference between T^ a v and 
G^ a v /%-kG is of second order, T^ a v on the right hand side of Eq.(j5H]) can be replaced by G^ a v . For 
the v = component, this leads to 

r(%„ = ^(-<r( 2 ^ Q/ ,>G(°)%+<r( 2 )%^>G(°)%). (56) 

By calculation, 



a" a' 

<r^ l0 >=-h* k h' kJ , 
< r (2)fc a >= \tt l (h? k , x -h jhk -h khj ), 



one obtains 



= ^5-^ W - (57) 



8vrGa 2 

In this form , r' '''^ is a coupling between the gravitational waves and the background spacetime of 
expanding universe. Thus, Eq . ()52|) is cast into the form 

^ + m( Pgw + Pgw ) = -^-^ h ' kl h k \ (58) 

which is the same as Eq. (T45|) with 



V gw =P gw + ^f^^. (59) 

By now, the origin of the pressure V gw is clear. In presence of matter, the equation of energy 
conservation can still be written in a form of continuity equation as Eq . ff45T) . where, for the pressure, one 



has to use V gw to take into account the interaction between GW and matter. Since the term ^ G ^ a -i h' kl h kl 
as the difference V gw — p gw is a typical viscosity term for RGW since it involves the time derivative h' M . 
The prefactor 24nGa i determines the direction and rate of energy interchange between RGW and matter. 
But, for RGW whose wavelengths are within the horizon of the expanding universe, this term is a small 
correction to p gw ~ g^Ga 1 ^ki^ k l < Since the rates T~L, H'/H are smaller than the frequencies of RGW of 
interest, and therefore \h kl \ <C \h' kl \. In short-wavelength limit, the correction term is recognized 

as being one order of 0(X/R) higher than p gw . 

Moreover, for the inflationary stage with the scale factor a(r) oc 1/r 

n 2 - n' = o, (60) 

one has V gw = p gw . For a general expansion stage with a(r) oc r a , one has 

n 2 -n' a + i 

—}C- = — (61) 
n 



which decreases quickly with the expansion, and V gw effectively approaches to p gw . This is true for 
the radiation-dominated stage (a = 1) and the matter-dominated stage (a = 2). Putting all these 
consideration together, for whole history of cosmic expansion since inflation, the difference V gw — p gw is 
negligibly small, and one can amply use Eq . fj50|) for energy conservation in practical studies of RGW. Only 
in the pre-inflationary stage the difference V gw —p gw can be significant for RGW of long wavelengths [TTJ. 

6. Effective EMPT and Wave Equation 

For the effective EMPT, we now explicitly show that the energy conservation is equivalent the field 
equation of RGW in the expanding RW spacetime. The field equation for RGW is 

h'la + 2Wi kl - d m d m h kl = 0, (62) 



which is commonly used in literature. The equation of energy conservation is given in Eq. (l45|) . Taking 
time derivative of p gw as given in Eq. PTj) . we have 

dp gw _ 2a' , I am h kl , u' h 'kh 



dr 8vrGa 3 



[Uh kl h kl + -(d m h kl d m h kl + h kl h kl )\ 



' [n'h' kl h kl + m' H h kl + Hh' kl h' kl + hd m h' kl d m h kl + huh'" 1 )} 



87rGa 2L kL Kl kl Kt ki kl 4 
-2Hp gw + 87rGa 2 W-'h' k ih k i + Hh kl h k i + Hh kl h' M 
1 



+-(d m h kl d m h kl + h kl h k % 



Using Eq.flUJ) for p gw and Eq. (l46]l for V gw to calculate %p gw + 2>l-LV gw yields 

H 2 — 7i' 

Tipgw + 'i'WPgw = Hpgw + ?>Hpg W H g 7r (^ cl 2 ^2 

1 ;[(2u 2 -n')h M h kl + ud m h M d m h kl 



(63) 



8vrGa 2 1 

~Hh' kl h' kl ]. (64) 

Adding these together, by using d m h' kl d m h kl = —h' kl d m d m h kl plus a surface term which is dropped, we 
finally arrive at the following result: 

^ + m( Pgw + V gw ) = ^2(\h' M + m kl )(h' kl + 2Hh' kl - d m d m h kl ). (65) 



From this equality, one sees that the wave equation (|62|) implies the equation of energy conservation Q45J) , 
and vice versus. Considering that for most of stages of cosmic expansion V gw ~ p gw , one can amply use 
either the wave equation (}6"2l) or the conservation equation (i50l) to study RGW. 
In the special case of % = 0, Eq. (i65l ) reduces to (choosing a = 1) 



dr 32ttG 
the same as Eq. (fT0j) for Landau-Lifshitz' case. 



h kt (h kl -d m d m h kl ), (66) 
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7. Conclusion 



The wave equation of GW is shown to be equivalent to the equation of energy conservation for 
Landau-lifshitz' definition of the pseudo-tensor t a @ in flat spacetime. It is explicitly demonstrated that, in 
the short wavelength limit after Brill-Hartle average, both Landau-lifshitz' t a/3 and the effective t"^ lead 
to the gauge-invariant, concise expression t^jg = < W ^h^p >, which is often used in study of 
gravitational waves. 

For RGW in the RW spacetime, the two forms of pressures p gw and V gw are analyzed, and their 
difference is traced to be V gw — p gw = g^T^ ^ | M - When matter is present, the energy conservation of 
RGW can still be written as a form of continuity equation + 3T-L(p gw + V gw ) = 0, where the pressure 
V gw takes into account the interaction between RGW and matter. For wavelengths within the horizon, 
the difference V gw —p gw is small, compared to p gw , and decreases with the cosmic expansion. During the 
de Sitter stage V gw - p gw = 0. 

We have also explicitly demonstrated that the wave equation ri' kl + 21-Ui kl — d m d m hu = is equivalent 
to %^ + 3H{p gw + V gw ) = 0. Nevertheless, since V gw ~ p gw after inflation, one can amply use either 
the wave equation or the "canonical" equation + 3H(p gw +p gw ) = in study of RGW for practical 
purposes. 
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